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Abstract 

The ground state properties of the two dimensional spatially anisotropic 

Heisenberg model are investigated by use of field theory mappings, spin- wave 

expansion and Lanczos technique. Evidence for a disorder transition induced 

by anisotropy at about Jy/Jx < 0.1 is shown. We argue that the disordered 

phase is gapless and its long wavelength properties can be interpreted in terms 

of decoupled one dimensional chains. 
75.10.-b, 75.30.Ds,75.10.Jm 
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The search for disordered (spin hquid) ground states in two dimensional (2D) electronic 
models has been pursued since the seminal work of Fazekas and Anderson [|T| on quantum 
antiferromagnets (QAF) in frustrating lattices. This problem has been revived in the last 
few years due to the resonating valence bond conjecture 0] which stimulated much numerical 
work on the subject. Despite the considerable effort, the existing evidence in favor of a spin 
liquid in 2D frustrated QAF is weak at best with the single exception of the Kagome 
QAF where a disordered ground state is plausible even if the presence of non conventional 
magnetic order is still a possibility. 

In this Letter, we present analytical as well as numerical evidence supporting an order- 
disorder transition in the square lattice S = 1/2 QAF driven by spatial anisotropy in 
the nearest neighbor coupling. This model does not introduce frustration and therefore 
presents several advantages with respect to the previously investigated systems, the most 
important being the absence of any plausible order parameter competing with the Neel 
staggered magnetization m = J2r Sr, exp(Q ■ R) [Q = (vr, vr)]. The model is defined by the 
hamiltonian 

(1) 

R 

where Sr are spin 1/2 operators living on a square lattice, x and y are unit vectors and 
^ Jy ^ Jx- The isotropic limit {Jy = J^) has been extensively studied by exact diagonal- 
izations and quantum Monte Carlo |^,|| with the resulting evidence of a finite staggered 
magnetization in the thermodynamic limit [|| m ~ 0.3075 quite close to the spin-wave theory 
(SWT) estimate m = 0.3034 0. Physically, the strongly anisotropic model (P describes a 
system of weakly coupled AF chains whose study has attracted considerable interest among 
theoreticians and experimentalists in view of the possibility to observe the peculiar features 
of one dimensional physics |]10| . 



The presence of an order-disorder transition in model ([1|) has been conjectured by several 
authors |n],|T2[ and can be motivated by the standard mapping of the 2D quantum model 
(0) into the (2+1) dimensional 0(3) non linear sigma model (NLaM) defined by the action 



S = ^J dxdydt [r^id^nf + Tyidyuf + Xo{dtnf] (2) 

where n is a unit vector. The lowest order estimates of the parameters give T^^ = Jx/4:, 
'^y — Jyl^i Xo^ = ^ci^iJ'x + Jy) where a is the lattice spacing. Two limits of the action 
(@) can be easily analyzed: The isotropic model is known to be ordered for the physically 
relevant parameters ||l3l while the Jy ^ Q limit of the action (H) correctly describes a stack 
of uncoupled (1+1) dimensional models which are disordered at any finite "temperature" 
g = (T^xoCf^)^^^^ owing to the Mermin Wagner theorem. Most interestingly, the order- 
disorder transition occurs at a finite value of the spatial anisotropy and belongs to the 
universality class of the classical three dimensional Heisenberg model. 

Similar conclusions can be drawn directly from SWT on the quantum hamiltonian (|l]). In 
fact, the 1/S expansion of the staggered magnetization can be straightforwardly generalized 
to anisotropic models and predicts a breakdown of Neel order at about ~ 0.03 (0.07) at 
first (second) order in 1/25*. The increase of the critical anisotropy parameter ac in going 
from first to second order gives confidence about the actual occurrence of the transition, 
which is in fact enhanced by quantum fluctuations. Therefore, on the basis of the field 
theory mapping and SWT, we expect that by lowering the anisotropy parameter a = Jy/ Jx 
a disordered phase sets in within a finite interval > a > 0. This prediction should be 
qualitatively correct because field theory methods are known to reproduce the physics of 



QAF both in the isotropic two dimensional limit ||T3[ and in the one dimensional {a = 0) 



case, provided the topological term is included in (^ |T^ . 

In order to test the theoretical predictions on the model (|l]) and to determine the prop- 
erties of the two dimensional spin liquid state, we have performed Lanczos diagonalizations 
on several lattices up to 32 sites. A proper finite size scaling of small lattice results is ob- 
viously important in order to provide a correct interpretation of the diagonalization data. 



Fortunately, much work has been recently devoted to this subject [|T3|,|TB[ showing that, in 



the ordered phase, the renormalization group flow drives the model towards weak coupling 
making SWT asymptotically exact at long wavelengths. As a consequence, SWT is able 
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to describe the leading size corrections in finite lattices and therefore represents a powerful 
method for analyzing small size data, particularly in non frustrated models. SWT, general- 
ized to finite systems [jl^ , compares quite favorably to numerical results in bipartite lattices. 
The good agreement also persists for the anisotropic model (|I|), as shown in Fig. 1 where a 
finite size estimate of the order parameter is plotted as a function of the anisotropy a. The 
breakdown of SWT at a < 0.1 again suggests that a qualitative change in the ground state 
properties is occurring in this regime. 

Further numerical evidence of the phase transition can be obtained by the structure 
of the energy spectrum as a function of the total uniform magnetization. According to 
a recent analysis the presence of Neel long range order in the thermodynamic limit 
reflects in the structure of the energy spectrum in finite size systems. In fact, if long range 
antiferromagnetic order is present in the system, the dependence of the energy E{S) on the 
total spin S must follow the approximate relation: 

E{S) ^ E{0) + S{S + 1)/{2In) (3) 

up to a maximum value Smax of the order of the square root of the number of sites A^. Eq. 
approximately reproduces the energy spectrum of a spin-S* rigid rotator with a momentum of 
inertia per site In /N corresponding to the uniform susceptibility x of the model. Notice that 
this criterion correctly reproduces the absence of antiferromagnetic long range order in one 
dimension where the energy spectrum scales as [|l|] E{S) ~ E(0) + S^/{2xN). The relation 
(0) is in fact quite well verified both for the isotropic system and for the anisotropic ones 
up to a ~ 0.1 where significant discrepancies appear (see Fig. 2a). The strong deviations 
from Eq. (0) which develop in the numerical data can be directly related to the asymptotic 
decoupling of the chains in the square lattice, leading to an approximate linear dependence 
E{S) ~ -E(O) + AS*. This anomaly does not seem to scale to zero in the thermodynamic 
limit but instead persists in all the lattices we have analyzed. From our finite size data 
an accurate estimate of the uniform susceptibility can be extracted by a quadratic fit of 
the energy spectrum E{S) at small but finite (uniform) magnetization S/N. The results 
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can then be extrapolated to infinite volume by use of a finite size scaling of the form |^ 
= + AN~^^'^ + B/N which is known to work both in the isotropic limit and in 
one dimension (extreme anisotropy). The extrapolation, shown in Fig. 2b for several values 
of the anisotropy parameter, is in fact in good agreement with the known values at a = 
and a = 1 (also shown in the figure) and predicts a smooth, featureless susceptibility in the 
whole anisotropy range. This result indicates that the system remains gapless across the 
phase transition and suggests that the nature of the disordered phase of this model might 
be more exotic than the expected nondegenerate singlet, in agreement with a conjecture put 
forward by Haldane ||14|| . 

However, the possibility of a gapless phase contrasts with the commonly accepted phase 
diagram of the model (|l|) defined on two chains |T8| , p!9[| where a gap A{a) is believed to open 



at every finite value of the anisotropy parameter a. The disordered phase in the two chain 
model is in fact continuously connected with the a — > oo limit where the gap is interpreted 
as the effect of the finite size of the lattice along the y direction. The same result actually 
holds for every even number of chains while an odd chain model remains gapless all the way 
to the a oo limit. Therefore it is not too surprising that the hamiltonian (|T]) on square 
clusters preserves the peculiarities of the odd chain sequence and does not open a gap at any 
a. In finite clusters, however, a gap is always present and we must investigate whether it 
disappears in the thermodynamic limit. We have analyzed the finite size scaling of the gap 
in the case of two chains (L x 2), three chains {L x 3) with antiperiodic boundary conditions 
along the y direction, and for square clusters. For any a the lowest excited state is always a 
triplet but its size dependence is quite different in the three cases. In order to see whether 
a gap is present in the strong anisotropic region we assumed that, for a — and L — > oo, 
the gap A(L, a) can be expressed in a scaling form, as usual near a critical point: 

A(L,a) = A(L,0)F [aL(log L)^/2] (4) 

where the one dimensional gap A(L,0) is known to scale as 1/L. The specific form (^ 
has been chosen in order to match with first order perturbation theory in a and does not 



depend on the number of chains of our lattice. However, the scahng function F{x) behaves 
quite differently in the three geometries, as can be seen in Fig. 3. The correctness of our 
scaling form (^) can be inferred by the collapse of the finite size numerical data on a smooth 
curve in all cases, provided a is sufficiently small. The region where the universal curve 
F{x) is defined increases with growing size and the thermodynamic limit at fixed (small) a 
corresponds to the large x region of the scaling curve which should be extrapolated from 
the finite size data. In the two chain model F{x) clearly goes through a minimum and then 
grows, suggesting a linear asymptotic behavior at large x which implies a finite gap of order 
Jy at small a in agreement with field theoretical analysis UTSf. Instead, the scaling function 
is always monotonic both in the three chain case and, even more convincingly, in the square 
clusters, supporting the absence of a gap in these systems. 

In order to understand how a disordered gapless phase may appear in 2D it is useful to 
consider other physical quantities like the spin-wave velocity and the momentum dependence 
of the magnetic structure factor. Again, SWT provides a valuable help in the interpretation 
of the numerical results. The spin velocity is almost constant at all anisotropics ranging 
between the one dimensional value Cx = vr/2 and the isotropic limit ~ 1.56 which 

are both reproduced within 10% by second order SWT generalized to anisotropic models. 
A surprising result of SWT is the enhancement of the anisotropy in the spin velocity ratio 
induced by quantum fiuctuations: 

J)' =«^(«) Z{a) = l-^C{a) (5) 

^ , , 1 x-^ , (cos kx — COS ky) (cos k^ + a cos ky) 
^ k \/ (1 + aV — (cos kx + a cos kyY 



In fact, while at lowest order the spin velocity ratio coincides with the anisotropy parameter, 
the one loop calculation always reduces the Z{a) factor. Obviously, the correction C{a) 
vanishes at the isotropic point a = 1 but diverges logarithmically in the a — > limit. 
Therefore, SWT suggests the occurrence of a decoupling transition at a finite value of a 
signaled by Z{ac) = 0. The same anisotropy renormalization factor Z{a) governs the long 

6 



wavelengths behavior of the physical correlation functions. In particular, the magnetic 
structure factor behaves as 



S{kx,ky) (X + aZ{a)k'^ (6) 

In order to verify these predictions we tested Eq.(^ against Lanczos diagonalizations in the 
32 sites lattice. The results are shown in Fig. 4 together with the zero and one loop SWT 
results for the spin velocity ratio in the thermodynamic limit. The numerical data are in 
good agreement with the spin-wave results in the 32 site lattice, and show an even larger ef- 
fect. Therefore we are led to conclude that at long wavelengths a decoupling transition may 
actually occur in strongly anisotropic spin models. The phase diagram of the anisotropic 
model (0) suggested by SWT is depicted in Fig. 5 for generic spin S systems. The transi- 
tion line where the staggered magnetization vanishes has been calculated at the lowest order 
spin-wave level together with the locus Z{a) = where we expect a "decoupling transi- 
tion". At the same order in 1/5 we have found that these two lines approximately coincide 
up to a critical value of a beyond which the system disorders without decoupling. In the 
strong anisotropy limit the transition is characterized by the vanishing of both the staggered 
magnetization and the spin velocity ratio leading to a picture of basically uncoupled chains 
with interesting experimental consequences about the possibility to observe ID behavior in 
real systems. We believe that this phase diagram is qualitatively correct although higher 
order terms in the SWT expansion (available only for the magnetization) may quantitatively 
change the phase transition line. In order to fully characterize the disordered phase, topo- 
logical defects must be taken into account leading to a possible difference between integer 



and half integer spin systems [14 
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FIGURES 

FIG. 1. Order parameter m = ^S{'k, tt)/V for different lattice sites compared with the first 
(dashed hnes) and second order (continuous hne) SWT results. The Lanczos data are obtained on 
tilted lattices L\/2 x L^/2 with L = 2 (triangles), L = 3 (squares) and L = 4 (circles). The lowest 
curves refers to the infinite size SWT results. 

FIG. 2. a): rigid rotator anomaly 5x~^ = y{E{S) - E{Q)) 

— 1/{2xl), where xl is estimated on finite sizes for 8 (triangles), 18 (squares) and 32 (circles) 
lattice sites. The continuous lines are guides to the eye, the dashed line represents the ideal be- 
haviour in a quantum antiferromagnet. b): Finite size scaling of the inverse susceptibility. The 
finite size data (notation as in Fig. 1) extrapolated (see text) to infinite size (continuous line). The 
stars are the exact values in the isotropic and one dimensional case. 

FIG. 3. Diagonalization data of the gap scaling function (see Eq. ^) for the two chain model, 
three chains and square clusters. For the two and three chains, open triangles refer to L = 4, full 
triangles L = 6, open squares L = 8, full squares L = 10, open circles L = 12. For the square 
clusters open triangles correspond to 8 sites, full triangles to 18 sites, open squares to 32 sites and 
stars to the 4x4 cluster. 

FIG. 4. Square of the spin velocity ratio vs. anisotropy. The dashed line is the leading SWT 
result, the continuous line includes the one loop correction in the thermodynamic limit. Finite size 
estimates on a 32-site lattice for aZ{a) are obtained by exact diagonalization (full triangles) and 
second order finite size SWT |jT^ (open triangles). 

FIG. 5. Phase diagram of the spatially anisotropic Heisenberg model obtained via one loop 
SWT. The order parameter vanishes along the continuous line and the spin-wave velocity ratio 
along the dashed line. The long dashed line indicates a crossover transition between a decoupled 
phase (DP) and a normal disordered phase with a finite spin velocity ratio. 
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